ABSTRACT
COVID-19 resulted in extensive morbidity and mortality worldwide. SARS-CoV-2 evolved rapidly, with increasing transmission due to Variants of Concern (VOC). Identifying VOC became important but genome submissions from low-middle income countries (LMIC) remained low leading to gaps in genomic epidemiology. We demonstrate the use of a specific mutation RT-PCR based approach to identify VOC in SARS-CoV-2 positive samples through the pandemic in Pakistan. We selected 2150 SARS-CoV-2 PCR positive respiratory specimens tested between April 2021 and February 2022, at the Aga Khan University Hospital Clinical Laboratories, Karachi, Pakistan. Commercially available RT-PCR assays were used as required for mutations in Spike protein (N501Y, A570D, E484K, K417N, L452R, P681R and deletion69_70) to identify Alpha, Beta, Gamma, Delta, and Omicron variants respectively. Three pandemic waves associated with Alpha, Delta and Omicron occurred during the study period. Of the samples screened, VOC were identified in 81.7% of cases comprising mainly; Delta (37.2%), Alpha (29.8%) and Omicron (17.1%) variants. During 2021, Alpha variants were predominant in April and May; Beta and Gamma variants emerged in May and peaked in June; the Delta variant peaked in July and remained predominant until November. Omicron (BA.1) emerged in December 2021 and remained predominant until February 2022. The CT values of Alpha, Beta, Gamma and Delta were all significantly higher than that of Omicron variants (p<0.0001). We observed VOC through the pandemic waves using spike mutation specific RT-PCR assays. We show the spike mutation specific RT-PCR assay is a rapid, low-cost and adaptable for the identification of VOC as an adjunct approach to NGS to effectively inform the public health response. Further, by associating the VOC with CT values of its diagnostic PCR we gain information regarding the viral load of samples and therefore the level of transmission and disease severity in the population.
ABSTRACT
BACKGROUND: We investigated the genome diversity of SARS-CoV-2 associated with the early COVID-19 period to investigate evolution of the virus in Pakistan. MATERIALS AND METHODS: We studied ninety SARS-CoV-2 strains isolated between March and October 2020. Whole genome sequences from our laboratory and available genomes were used to investigate phylogeny, genetic variantion and mutation rates of SARS-CoV-2 strains in Pakistan. Site specific entropy analysis compared mutation rates between strains isolated before and after June 2020. RESULTS: In March, strains belonging to L, S, V and GH clades were observed but by October, only L and GH strains were present. The highest diversity of clades was present in Sindh and Islamabad Capital Territory and the least in Punjab province. Initial introductions of SARS-CoV-2 GH (B.1.255, B.1) and S (A) clades were associated with overseas travelers. Additionally, GH (B.1.255, B.1, B.1.160, B.1.36), L (B, B.6, B.4), V (B.4) and S (A) clades were transmitted locally. SARS-CoV-2 genomes clustered with global strains except for ten which matched Pakistani isolates. RNA substitution rates were estimated at 5.86 x10-4. The most frequent mutations were 5' UTR 241C > T, Spike glycoprotein D614G, RNA dependent RNA polymerase (RdRp) P4715L and Orf3a Q57H. Strains up until June 2020 exhibited an overall higher mean and site-specific entropy as compared with sequences after June. Relative entropy was higher across GH as compared with GR and L clades. More sites were under selection pressure in GH strains but this was not significant for any particular site. CONCLUSIONS: The higher entropy and diversity observed in early pandemic as compared with later strains suggests increasing stability of the genomes in subsequent COVID-19 waves. This would likely lead to the selection of site-specific changes that are advantageous to the virus, as has been currently observed through the pandemic.